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Abstract 

We study the impact of the symmetry energy on properties of nuclear pasta phases and crust- 
core transition in neutron stars. We perform a self-consistent Thomas-Fermi calculation employing 
the relativistic mean-field model. The properties of pasta phases presented in the inner crust of 
neutron stars are investigated and the crust-core transition is examined. It is found that the slope 
of the symmetry energy plays an important role in determining the pasta phase structure and 
the crust-core transition. The correlation between the symmetry energy slope and the crust-core 
transition density obtained in the Thomas-Fermi approximation is consistent with that predicted 
by the liquid-drop model. 
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I. INTRODUCTION 


Neutron stars are great laboratories for the study of asymmetric nuclear matter over 
a wide range of density. With increasing depth in the neutron star, the matter density 
may rise from the subnuclear region to several times normal nuclear matter density [l|, 
0]. In general, a neutron star consists of an outer crust of nuclei in a gas of electrons, 
nuclei in a gas of neutrons and electrons, and a liquid core of uniform 
4], The inner crust of neutron stars has drawn much attention due to 


an inner crust o 


dense matter 


its complex phase structure and significant role in astrophysical observations (sUsl. In the 
inner crust, spherical nuclei may become unstable as the density increases toward the crust- 
core transition, and the stable nuclear shape is likely to change from droplet to rod, slab, 


tube, and bubble, known as nuclear pasta phases 


10]. The crust-core transition occurs 


at the density where the energy density of the homogeneous matter becomes lower than 
that of the pasta phase. During the last decades, the properties of pastaphases have been 
investigated by using various methods, such as the liquid-drop model [g, |l0|, [ll| and the 
Thomas-Fermi approximation jg, 9, 12, 131. In Ref. 12], the density region of nonspherical 
nuclei was evaluated by using a parametrized Thomas-Fermi approximation, which was 
found to be sensitive to the density dependence of the nuclear symmetry energy. In Ref. [;8], 
a self-consistent Thomas-Fermi approximation was used to calculate properties of the inner 
crust including pasta phases, and it was found that the symmetry energy and its slope 
could have significant impacts on the pasta phase structure and crust-core transition. In 
our previous work 10], the effects of the symmetry energy on pasta phase properties and 
crust-core transition were investigated by employing the coexisting phases method based 
on a liquid-drop model, and the correlation between the symmetry energy slope and the 
crust-core transition was obtained and analyzed. It is noticeable that the symmetry energy 
and its slope could play an important role in determining the pasta phase structure and 
crust-core transition in neutron stars. 

In recent years, the nuclear symmetry energy and its density dependence have received 


great interest due to 
and astrophysics 


14 


reir importance for understanding many phenomena in nuclear physics 


17]. The value of the symmetry energy E sym at saturation density 


is constrained by experiments to be about 30 ± 4 MeV, while its slope L at saturation 
density is still very uncertain and may vary from about 20 to 115 MeV jisl ]. It has been 
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found that various properties of neutron stars, such as the crust structure, the crust-core 
transition, and the star radius, are sensitive to the symmetry energy E sym and its slope 


L [d, 12, 19, 3- In Ref. 


19], the impact of the symmetry energy on the crust-core transition 


was examined using various effective Skyrrne and relativistic approaches, in which the crust- 
core transition density obtained from the dynamical and thermodynamical methods showed 
a clear decrease with increasing L. The correlation between the crust-core transition density 
and the symmetry e nerg y slope L has been extensively studied in the literature using various 
methods |8|, [l0|, [12, |19j]. It is shown that the resulting transition density depends on the 
method and effective nuclear interaction used in the calculation. It is important to make 
further investigations in order to determine a clear correlation between the symmetry energy 
slope L and the crust-core transition. 

The main purpose of this article is to investigate the impact of the symmetry energy 
on pasta phase properties and explore the correlation between the symmetry energy slope 
L and the crust-core transition. We perform a self-consistent Thomas-Fermi calculation 
employing the relativistic mean-field (RMF) model j^l, 22] for nuclear interactions. In 
the Thomas-Fermi approximation, the surface effect and nucleon distributions are treated 
self-consistently, rather than by assuming a sharp interface as in the coexisting phases 
method [lo|. In our most recent study [23], we made a detailed comparison between the 
Thomas-Fermi (TF) approximation and the coexisting phases (CP) method with only the 
droplet configuration. It is interesting to compare their difference in pasta phases includ¬ 
ing all configurations as mentioned above. For the nuclear interaction, we adopt the RMF 
model with two different parametrizations, TM1 24] and IUFSU [25|, both of which are 
known to be successful in describing the ground-state properties of finite nuclei, including 
unstable ones. In the RMF approach, nucleons interact via the exchange of scalar and 
vector mesons, and the model parameters are generally fitted to nuclear matter saturation 
properties or ground-state properties of finite nuclei. The TM1 parametrization includes 
nonlinear terms for both a and uj mesons, while an additional co-p coupling term is added in 
the IUFSU parametrization. It is well known that the u>-p coupling term plays a crucial role 
in modifying the density dependence of the symmetry energy and affecting the neutron star 


properties [16], 20, |25h27]. In order to evaluate the impact of the symmetry energy slope L 


on pasta phase properties and crust-core transition, we employ two sets of generated models 
based on the TM1 and IUFSU parametrizations as given in Ref. 23]. The model parameters 
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were determined by simultaneously adjusting g p and A v so as to achieve a given L at satu¬ 
ration density % while keeping E sym fixed at a density of 0.11 frn -3 . We note that the fixed 
density in Ref. 10] was chosen to be the saturation density, namely, ng x = n 0 , whereas ng x 


23] - It has been shown in Ref. 23] that the choice of the fixed 


= 0.11 fm -3 was used in Ref. 
density ng x = 0.11 fm -3 could produce very similar binding energies for finite nuclei within 
one set of generated models. Furthermore, all models in each set have the same isoscalar 
saturation properties and fixed symmetry energy at Ufi x = 0.11 fm -3 , but they have different 
symmetry energy slope L. By using the set of models with different L, it is possible to study 
the impact of L on pasta phase properties and explore the correlation between L and the 
crust-core transition. 

This article is organized as follows. In Sec. HH we briefly describe the RMF model and the 
self-consistent TF approximation used in this study. In Sec. IIII1 we present the numerical 
results and examine the impact of the symmetry energy on pasta phase properties, while the 
correlation between the symmetry energy slope L and the crust-core transition is discussed. 
Section m is devoted to the conclusions. 

II. FORMALISM 


7 approximation by employing the 


,[22], nucleons interact through 


The inner crust of neutron stars is studied within the T. 

RMF model for nuclear interactions. In the RMF model |21 
the exchange of various mesons. The mesons considered here are the isoscalar-scalar meson 
er, the isoscalar-vector meson u, and the isovector-vector meson p. For a system consisting 
of protons, neutrons, and electrons, the Lagrangian density reads 


C 


RMF 
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( 1 ) 


where W^ u , R afiu , and F^ are the antisymmetric field tensors corresponding to u; M , 
and respectively. In the RMF approach, the meson fields are treated as classical fields, 
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and the field operators are replaced by their expectation values. For a static system, the 
nonvanishing expectation values are a = (a), u ; = (cj°), p = (p 30 ), and A = (A 0 ). The 
equations of motion for these mean fields derived from the Lagrangian density dU) have the 


following form: 

-W + mla + g 2 a 2 + g 3 a 3 = -g a ( n p + <) , (2) 

-V 2 (n + mluj + c 3 u> 3 + 2A v p 2 p 2 p 2 (n = g u (n p + n n ), (3) 

- V 2 p + m 2 p p + 2A v glg 2 p uj 2 p = y (n p - n n ), (4) 

—V 2 A = e (n p - n e ), (5) 


where nf and n, denote, respectively, the scalar and number densities of species i. The 
equations of motion for nucleons give the standard relations between the densities and 
chemical potentials, 


p,p — yJ (kip) + M* 2 + g u u: + y p + eA, 
p n = \J ( k n F f + M* 2 + g u u: - ^p, 


( 6 ) 

( 7 ) 


where M* = M + g a a is the effective nucleon mass, and k l F is the Fermi momentum of 
species i, which is related to the number density by n t = (k F ) 3 /37r 2 . 

The matter in the inner crust of neutron stars contains protons, neutrons, and electrons 
under the conditions of (3 equilibrium and charge neutrality. We employ the Wigner-Seitz 
cell approximation to describe the inner crust, in which the equilibrium state is determined 
by minimization of the total energy density at zero temperature. The stable cell shape may 
change from droplet to rod, slab, tube, and bubble as the density increases. For simplicity, we 
assume the electron density is uniform throughout the Wigner-Seitz cell, since the electron 
screening effect is known to be negligible at subnuclear densities 28]. Furthermore, we also 


neglect the correction caused by the Coulomb interaction with charged particles in other 


cells, which is negligibly small in most cases [29 
energy per cell is calculated from 


30]. In the TF approximation, the total 


-Fceii / ^nnr(r)dr T E e VAn, 


( 8 ) 


' cell 


where e e is the electron kinetic energy density, and £ rm f(r) is the local energy density at 
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position r, which is given in the RMF model by 

r -k 


£rmf(r) = “2 f k 2 Vk 2 + M ’ 
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+i(V<x ) 2 + \m 2 a o 2 + \g 2 cr 3 + 
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n P 


(9) 


Here, we consider different pasta configurations, including the droplet, rod, slab, tube, and 
bubble. The volume of the Wigner-Seitz cell for different configurations can be written as 


Kell — < 


4 S 

3 n 1 ws 

(droplet and bubble), 

/yrr 2 

1 ws 

(rod and tube), 

2 r / 2 

ws 1 

(slab), 


( 10 ) 


where r ws is the radius of a spherical cell for the droplet and bubble configurations, while 
the rod and tube have cylindrical shapes with radius r ws and length l , and the slab has 
width l and thickness 2r ws . We note that the choices of the length for a cylindrical shape 
and the width for a slab are somewhat arbitrary jg], which would not affect the resulting 
energy density of the system. 

At a given average baryon density n&, we minimize the total energy density with respect 
to the cell size r ws for each pasta configuration, and then we compare the energy densities 
between different configurations in order to determine the most stable shape that has the 
lowest energy density. Furthermore, the energy density of the corresponding homogeneous 
phase at the same 7i b is also computed, and the crust-core transition occurs at the density 
where the energy density of the homogeneous phase becomes lower than that of the pasta 
phase. In order to calculate the total energy per cell given by Eq. (jHJ) at fixed r ws and n b , 
we solve the coupled Eqs. (J2J) (EJ) under the constraints of (3 equilibrium, charge neutrality, 
and baryon number conservation, which have the following form: 


hn hp T [7ei 


N e = N p = 


n p (r)dr, 


'cell 


nbVceii = / [n p (r) + n„(r)]dr. 


( 11 ) 

( 12 ) 

(13) 


' cell 
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In practice, we start with an initial guess for the mean fields cr(r), (u(r), p{ r), and A(r), then 
determine the chemical potentials /i n , p p , and p e by the constraints (fTT]) -(fT3 |) . Once the 
chemical potentials are obtained, it is easy to calculate various densities and solve Eqs. (I?l) - 
([5]) to get new mean fields. This procedure should be iterated until convergence is achieved. 


III. RESULTS AND DISCUSSION 


In this section, we present numerical results for the inner crust of neutron stars, and 
we discuss the impact of the symmetry energy on pasta phase properties and crust-core 
transition. The results obtained from the self-consistent TF calculation are compared with 


those obtained using the CP method 
two successful RMF models, TM1 [24 


3 - 


For the effective nuclear interaction, we consider 


and IUFSU 


25]. The parameter sets and saturation 


properties of these two models are given in Tables HI and UTl respectively. In order to clarify the 
correlation between the symmetry energy slope L and the crust-core transition, we employ 
two sets of generated models based on the TM1 and IUFSU parametrizations as given in 
Ref. 23]. It is noticeable that all models in each set have the same isoscalar saturation 
properties and fixed symmetry energy E sym at a density of 0.11 fnU 3 but have different 


symmetry energy slope L. These models have been generated by simultaneously adjusting 
g p and A v so as to achieve a given L at saturation density n 0 while keeping E sym fixed at a 


density of 0.11 fm -3 as described in Ref. 23]. The parameters, g p and A v , generated from 
the TM1 and IUFSU models for different L are given in Tables IIIII and IIV1 respectively. In 
Fig. [Q we plot the symmetry energy E syra as a function of the baryon density rib for the 
two sets of models generated from TM1 (upper panel) and IUFSU (lower panel). One can 
see that all models in each set have the same E sym at a density of 0.11 fm -3 , but they have 
different values of E sym at lower and higher densities due to the difference in the slope L. It 
is obvious that a smaller L corresponds to a larger (smaller) E sym at lower (higher) densities. 
It will be shown below that the behavior of E sym plays a crucial role in determining the pasta 
phase structure and the crust-core transition. 

We first present the phase diagram for the inner crust of neutron stars and discuss the 
influence of the symmetry energy slope L on the pasta phase structure. In Fig. [2J, the density 
ranges of various pasta phases obtained from the self-consistent TF calculation are displayed 
for the two sets of generated models, IUFSU (left panel) and TM1 (right panel). It is found 
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that only the droplet configuration can occur before the crust-core transition for L > 80 
MeV, whereas the pasta phase structure may change from droplet to rod, slab, tube, and 
bubble for smaller values of L (e.g., L = 50 MeV). As one can see from Fig. [21 the onset 
density of nonspherical nuclei, i.e., the transition density from droplet to rod, significantly 
decreases with decreasing L in the low-L region (L < 70 MeV). This behavior is consistent 
with that reported in Ref. 12]. It can be understood from a fission-like instability of spherical 


nuclei predicted in the liquid-drop model 


12 


31], in which a spherical liquid drop becomes 


unstable to quadrupolar deformations when the volume fraction of the liquid drop reaches 
the fission-instability criterion, u = ( r d /r ws ) 3 = 1/8, with r d and r ws being, respectively, 
the radii of the droplet and the Wigner Seitz cell. From Fig. [I] below, we can see that a 
smaller L corresponds to a larger value of r d /r ws in the droplet phase at low densities, and, 
therefore, the model with a smaller L results in an earlier onset of nonspherical nuclei, as 
shown in Fig. [21 On the other hand, the onset density of homogeneous matter significantly 
decreases with increasing L , which has also been observed in earlier studies js, 10, 12, 17 ]. 
This trend may be understood from the energy-density curvature of pure neutron matter. 
According to the analysis in the liquid-drop model 17], the energy-density curvature of 
pure neutron matter at saturation density, Onm(^o)? is approximately proportional to L, 
and the inhomogeneous phase occurs in the spinodal region of asymmetric nuclear matter 
where the energy density has a negative curvature. A rough analysis implies that the larger 
Cnm(^o) is, the farther away a spinodal border of ^-equilibrium matter deviates from the 
saturation density no. Therefore, the model with a larger L leads to a smaller onset density of 
homogeneous matter, as shown in Fig. [2] [see also Fig.[9](a)]. We conclude that a smaller value 
of L results in an earlier onset of nonspherical nuclei and a later transition to homogeneous 
matter. Hence, the model with a smaller L predicts a more complex phase structure and a 
larger density range of pasta phases in neutron star crusts. 

The phase diagram obtained in the present TF calculation is very similar to that obtained 
using the CP method ID]. However, the bubble configuration could not appear even with 
the lowest L by using the CP method [hJ. In order to make a detailed comparison of the 
pasta phase structure between the TF and CP calculations, we present in Table M the onset 
densities of pasta phases for the original IUFSU (L = 47.2 MeV) and TM1 (L = 110.8 MeV) 
models. It is shown that only the droplet configuration appears in the case of TM1 for both 
TF and CP, which is due to its very large value of L. On the other hand, the original IUFSU 
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model predicts that almost all pasta phases can occur before the crust-core transition due to 
its small value of L. From Table IVl we can see that there are visible differences in the onset 
densities between TF and CP. This is mainly caused by the different treatments of surface 
and Coulomb energies 


238 ], which play an essential role in determining the phase shape. It is 
found that the onset density of homogeneous matter, i.e., the crust-core transition density, 
obtained in the TF approximation is slightly higher than that obtained using the CP method. 
This is because the configuration space of TF is much larger than that of CP, and, therefore, 
a lower energy density can be achieved in the minimization procedure of the TF calculation, 
which leads to a higher crust-core transition density. In Fig. [3] we plot the energy per 
nucleon of the pasta phase relative to that of homogeneous matter, A E, as a function of 
the average baryon density, rib, obtained from the TF and CP calculations using the original 
IUFSU model. One can see that there are significant differences in A E between TF and CP 
at lower densities, whereas the differences become much smaller in the pasta phase region. 
The comparison between TF and CP at low density with the droplet configuration has been 
extensively discussed in our previous work 23], and it was found that the simple CP method 
could not describe the nonuniform matter around the neutron drip density due to its energy 
being higher than that of homogeneous matter. From Fig. [3] it is seen that A E of the CP 
calculation becomes positive at rib < 0.003 frn -3 , which is consistent with the onset density 
of the droplet phase given in Table |V] On the other hand, the crust-core transition occurs 
at relatively high density, beyond which A E becomes positive. Although the differences in 
A E between TF and CP calculations at higher densities are rather small (on the order of 
a few keV), it may lead to visible differences in the onset densities as shown in Table [V] 
We note that the differences in the energy per nucleon between different pasta shapes are 
of the order of 0.1-1 keV, which has also been reported in Refs. 11, |29j. Therefore, the 
pasta phase structure is very sensitive to the method and nuclear interaction used in the 
calculation. 

In the present study, we focus on the correlation between the symmetry energy slope L 
and pasta phase properties. In Fig. [4] the size of the Wigner-Seitz cell, r ws , and that of 
the inner part (nucleus or hole), r d , which are obtained from the TF calculation with two 
extreme values of L in the TM1 (upper panel) and IUFSU (lower panel) sets, are displayed 
as a function of rib■ In the TF approximation, there is no distinct boundary between the 
liquid phase and the gas phase, so we prefer to define the size of the inner part, r d , by density 
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fluctuations, similar to that of Refs 

( n p) 

| 1 WS \ 

r d = 


.a 


281 ]. as 


' ws V H) 


( n p ) 2 


1 ID 


(droplet, rod, and slab), 
(tube and bubble), 


( 14 ) 


with D = 1,2,3 being the geometrical dimension of the system. The average values, (n p ) 
and (rip), are calculated over the cell volume V ce n given by Eq. (ITO]) . From Fig. [4j one can 
see that the model with the smallest L yields a rather complex structure of the inner crust, 
whereas only the droplet phase appears in the case with the largest L. It is found that r ws 
and show significant jumps at the transition points between different pasta shapes, which 
indicates that the transition is first order. One can see that r ws decreases with rib at lower 
densities, while it rapidly increases before the crust-core transition in the case of small L. 
A similar behavior of r ws was also observed in Ref. js]. 

It is of interest to compare the results between the TF and CP methods. In Fig. [51(a), 
r ws and r d obtained in the TF approximation are compared with those obtained by the 
CP method using the original IUFSU model. It is found that both r ws and from the 
CP method are somewhat smaller than those from the TF method. This can be explained 
by the behaviors of the electron chemical potential /i e shown in Fig. |5](b) and the average 
proton fraction Y p shown in Fig. [5](c). As discussed in Ref. 23], the CP method could yield 
larger /i e and Y p in comparison to those of TF [see Figs. [5](b) and[5](c)]. This is because the 
surface and Coulomb energies are treated perturbatively in the CP method, whereas they 
are included self-consistently in the minimization of the TF method. Due to larger /r e and 


Y p , the CP method would give rise to smaller r ws as s’ 


rown in Fig. [5](a), which could be 


23]. From Figs. [5](b) and[5](c), one 


explained by their correlation in the liquid-drop model 
can see that there are small jumps at rib ~ 0.079 fm -3 in the TF case, which are caused by 
the transition from slab to tube. It is obvious that trends in the differences between CP and 
TF are very similar for all pasta phases, and they are consistent with those obtained in the 
droplet phase j^fj]. 

To examine the influence of L on properties of the inner crust, we perform self-consistent 
TF calculations using the two sets of generated models based on the TM1 and IUFSU 
parametrizations. In Fig. [6] we present the following quantities: (a) the average proton 
fraction Y p , (b) neutron densities of the liquid phase and the gas phase, n nt L and n Uj G, at 
the center or boundary of the cell; and (c) the proton density of the liquid phase, n Pt L, at 
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the center or boundary of the cell. As one can see from Fig. [6](a), Y p decreases rapidly with 
increasing rib at lower densities and shows a significant L dependence at higher densities. It 
is found that a smaller L corresponds to a larger Y p at a fixed rib■ This trend is related to the 
density dependence of E sym shown in Fig. Q] Since E sym has been fixed at rib — 0.11 fm -3 , 


a smaller L in one set of generated models corresponds to a larger E sym at lower densities 
(rib < 0.11 fm -3 ). It is well known that a larger E sym favors a higher Y p in homogeneous 
/3-equilibrium matter. Therefore, a smaller L results in a larger Y p at the density close to 
the transition to homogeneous matter. The correlation between L and Y p is consistent with 
those reported in Refs, jg, 12]. A clear L dependence is also observed in Fig.^b), in which 
a smaller L corresponds to larger n n ^ and smaller n n) c■ The L dependence of n n L and 
n n ^G can be explained by the density dependence of E sym as discussed above. It is shown in 
Fig. [1] that a smaller L corresponds to larger E sym at rib < 0.11 fm -3 and smaller E sym at 
rib > 0.11 fm -3 . Therefore, the model with a smaller L favors more neutrons in the liquid 
phase and fewer neutrons in the gas phase, which results in larger n n> £ and smaller n n> G, as 
shown in Fig. [6])b). We note that the behaviors of n n j, and n n) c obtained in the present 


study are consistent with those reported in Refs. [8j, 12]. On the other hand, the behavior 


of n p i is somewhat complicated, as shown in Fig. [6])c). The model with a smaller L results 
in a more rapid decrease of n Pj L at 0.01 < rib < 0.05 fm -3 . This trend may be related to the 
LO-p coupling term and behaviors of chemical potentials as discussed in Ref. 23]. In Fig. 0 
we display chemical potentials of electrons, /x e , neutrons, p n , and protons, p p , as a function 
of rib obtained in the TF approximation for the two sets of generated models with several 
values of L. The chemical potentials are intensive quantities that play an important role in 
the TF calculation. With increasing rib, both p e and p n increase monotonically in all cases 
of L, whereas p p decreases. For the L dependence of chemical potentials at a fixed rib, it 
is found that the model with a smaller L generally results in larger p e and p n , as well as 
smaller p p . This can be explained by the contribution from the p meson in the RMF model. 
According to Eqs. ([6]) and (ITI) , the value of g p p plays an essential role in determining the 
chemical potentials p p , p n , and p e = //„ — p T> . One can see from Tables [HD and HV] that 
the model with a smaller L has relatively larger g p and A v , which yields a larger value of 
g p p. Therefore, the model with a smaller L leads to larger p e and p n , as well as smaller 
p p . By comparing results between TM1 and IUFSU, we find that the two sets of generated 
models have similar L dependence for all properties mentioned above. These results are also 
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consistent with those reported in Refs. [8|, |l0, 12, 23]. 

It is interesting to see how the density profile evolves through various pasta phases as 
the density increases. In Fig. [SI the density distributions of neutrons and protons inside 
the Wigner-Seitz cell are plotted at several values of rib- The calculations are performed 
with two extreme values of L in the set of IUFSU (L = 47.2 MeV and L = 110 MeV). As 
shown in Fig. [2J the original IUFSU model (L = 47.2 MeV) predicts that all pasta phases 
would occur before the crust-core transition, whereas only the droplet configuration appears 
in the case of L = 110 MeV. In Fig. [8] from top to bottom, we show the results of droplet, 
rod, slab, tube, and bubble, respectively. In the top panel at rib = 0.04 fm -3 , the results 
for a droplet with L = 47.2 MeV (thick lines) are compared to those with L — 110 MeV 
(thin lines). It is seen that the droplet size obtained with L = 47.2 MeV is larger than 
that with L — 110 MeV, and the nucleon distributions in the case of L = 110 MeV are 
more diffuse than those of L = 47.2 MeV. Moreover, as compared to the case of L = 110 
MeV, the original IUFSU model with L = 47.2 MeV yields lower neutron gas density at 
the boundary (n ni c) and higher neutron density at the center of the cell ( n nt L ), which are 
consistent with those shown in Fig. [6](b). In the second panel at rib = 0.06 frn -3 , the original 
IUFSU model predicts a rod phase, while the model with L — 110 MeV has undergone 


the crust-core transition. The transition from slab to tube occurs at rib ~ 0.079 fm 


-3 


m 


the original IUFSU model, and, therefore, sudden changes in the density distributions are 
observed by comparing the third and fourth panels; these cause visible jumps in chemical 
potentials and Y p as shown above. As one can see, from top to bottom, the distributions 
of neutrons and protons become more diffuse with increasing rib, and the neutron density 
difference between the liquid phase and the gas phase becomes smaller and smaller. The 
crust-core transition occurs at rib ~ 0.092 fm -3 in the original IUFSU model. 

Finally, we discuss the correlation between the symmetry energy slope L and the crust- 
core transition. In Fig. |9l we display the crust-core transition density rib,t and the proton 
fraction and the pressure at the transition point, Y p>t and Pt, as a function of L obtained 
from the self-consistent TF calculation using the two sets of generated models based on 
the TM1 and IUFSU parametrizations. As one can see from Fig. [9](a), there is a clear 


correlation between L and n^t, namely, rib,t decreases monotonica 


correlation is consistent with those reported in Refs. 8, h], 12, 19]. Compared to the results 
obtained using the CP method [loj], the transition densities obtained from the present TF 


ly with increasing L. This 
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calculation are slightly higher, which may be related to the differences in configuration space 
and treatments of surface and Coulomb energies. From Fig.[9](b), it is seen that Y P}t decreases 


significantly with increasing L. This trend is very similar to that observed in Refs. 10, [19]. 
It can be understood from the density dependence of E sym shown in Fig. [H As mentioned 
above, a smaller L in one set of generated models corresponds to a larger E sym at n b < 0.11 
fm -3 , and a larger E sym favors a higher Y p in homogeneous /3-equilibrium matter. Therefore, 
a smaller L results in a larger Y Pjt as illustrated in Fig. [9][b). The L dependence of the 
transition pressure Pt is nonmonotonic, as shown in Fig. |9][c). It is found that Pt decreases 
with increasing L in the large-L region (L > 60 MeV), whereas the opposite behavior is 
observed for L < 60 MeV. A similar behavior was also observed in Ref. [10 [.The nontrivial 

17]. For neutron- 


L dependence of P t is due to a competing effect, as discussed in Refs. [FOi, 



rich matter at fixed density and proton fraction, the pressure would increase with increasing 
L. However, the decrease of n b)t shown in Fig. |9][a) leads to a decrease of the pressure with 
increasing L. As a result, P t depends on L nonmonotonically, as illustrated in Fig. |9][c). 


IV. CONCLUSIONS 

In this study, we have investigated the impact of the symmetry energy on pasta phase 
properties within the self-consistent TF approximation. It has been found that the symmetry 
energy slope L plays an important role in determining the pasta phase structure and the 
crust-core transition. In order to clarify the influence of L, we have employed two sets 
of generated models based on the TM1 and IUFSU parametrizations that have the same 
isoscalar saturation properties and fixed symmetry energy at the density n b = 0.11 fm -3 but 
have different symmetry energy slope L. It has been observed that the model with a smaller 
L predicts an earlier onset of nonspherical nuclei and a later transition to homogeneous 
matter. For example, the original IUFSU model with L = 47.2 MeV predicts that the 
transition from droplet to rod occurs at n b ~ 0.049 fm -3 , then the pasta phases of slab, tube, 
and bubble appear one by one, and finally the transition to homogeneous matter occurs at 
n b ~ 0.092 fnu 3 . In contrast, the original TM1 model with L = 110.8 MeV predicts that 
only the droplet configuration appears at low density, and the transition from droplet to 
homogeneous matter occurs at n b ~ 0.062 fm -3 . Therefore, it can be concluded that the 
model with a smaller L results in a more complex phase structure and a larger density range 
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of pasta phases for neutron star crusts. In addition, some properties of the inner crust 
such as the proton fraction, chemical potentials, and density profiles have also been found 
to be correlated to the symmetry energy slope L , which could be partly explained by the 
density dependence of the symmetry energy. We have compared the results of the present 
TF calculation with those obtained using the CP method. It has been found that, although 
there are quantitative differences between these two methods, the qualitative behaviors of 
the inner crust are very similar to each other, and they are consistent with those reported 
in earlier studies BQ. 

The correlation between the symmetry energy slope L and the crust-core transition has 
been examined in the TF approximation using the two sets of models generated from the 
TM1 and IUFSU parametrizations. It has been found that the crust-core transition density 
rib,t decreases monotonically with increasing L. This correlation is consistent with those re¬ 
ported in Refs, js, 10, 12, 191. Also, the proton fraction at the transition point Y p t decreases 
significantly with increasing L, whereas the transition pressure Pt shows a nontrivial de¬ 
pendence on L. These correlations obtained from the present TF calculation are consistent 
with those obtained using the CP method 10]. In the TF approximation, nuclear shell and 
pairing effects have been neglected. It is of interest to further consider neutron paring and 
superfluidity in the inner crust of neutron stars. 
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TABLE I: Parameter sets used in this work. The masses are given in MeV. 


Model 

M m c m p 

9a 

9uj 

S P 

to 

? 

93 c 3 

A v 

TM1 

938.0 511.198 783.0 770.0 10.0289 

12.6139 

9.2644 

-7.2325 

0.6183 71.3075 

0.000 

IUFSU 939.0 491.500 782.5 763.0 

9.9713 

13.0321 

13.5900 

-8.4929 

0.4877 144.2195 0.046 


TABLE II: Saturation properties of nuclear matter for the TM1 and IUFSU models. The quan¬ 
tities Eq, K, E sym , and L are, respectively, the energy per nucleon, incompressibility coefficient, 
symmetry energy, and symmetry energy slope at saturation density no- 


Model n 0 (fm“ 3 ) E 0 (MeV) K (MeV) E sym (MeV) L (MeV) 


TM1 

0.145 

-16.3 

281.0 

36.9 

110.8 

IUFSU 

0.155 

-16.4 

231.0 

31.3 

47.2 


TABLE III: Parameters, g p and A v , generated from the TM1 model for different slope L at satu¬ 
ration density no with fixed symmetry energy E sym = 28.05 MeV at a density of 0.11 fm~ 3 . The 
last line shows the symmetry energy at saturation density, E sym (no). The original TM1 model has 
L = 110.8 MeV. 

L (MeV) 40.0 50.0 60.0 70.0 80.0 90.0 100.0 110.8 

gp 13.9714 12.2413 11.2610 10.6142 10.1484 9.7933 9.5114 9.2644 

A v 0.0429 0.0327 0.0248 0.0182 0.0128 0.0080 0.0039 0.0000 

E sym (n 0 ) (MeV) 31.38 32.39 33.29 34.11 34.86 35.56 36.22 36.89 
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TABLE IV: Parameters, g p and A v , generated from the IUFSU model for different slope L at 
saturation density no with fixed symmetry energy E sym = 26.78 MeV at a density of 0.11 fm“ 3 . 
The last line shows the symmetry energy at saturation density, E sym (no). The original IUFSU 
model has L = 47.2 MeV. 


L (MeV) 

47.2 

50.0 

60.0 

70.0 

80.0 90.0 100.0 110.0 

3p 

13.5900 

12.8202 

11.1893 

10.3150 9.7537 9.3559 9.0558 8.8192 

A v 

0.0460 

0.0420 

0.0305 

0.0220 

0.0153 0.0098 0.0051 0.0011 

E sy m(n 0 ) (MeV) 

31.30 

31.68 

32.89 

33.94 

34.88 35.74 36.53 37.27 


TABLE V: Comparison of the onset densities of pasta phases and homogeneous matter (horn.) 
between the TF and CP methods for the original IUFSU and TM1 models. 


Model 

Method 


Onset density (fm 3 ) 

droplet 

rod slab tube 

bubble hom. 

IUFSU 

TF 

0.000 

0.049 0.063 0.079 

0.085 0.092 

IUFSU 

CP 

0.003 

0.055 0.070 0.088 

0.089 

TM1 

TF 

0.000 

— — — 

0.062 

TM1 

CP 

0.001 

— — — 

0.058 
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FIG. 1: (Color online) Symmetry energy E sym as a function of the baryon density n^ for modified 
versions of TM1 (upper panel) and IUFSU (lower panel) with several values of L at saturation 
density. The symmetry energy is fixed at a density of 0.11 fm~ 3 . 
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FIG. 2: (Color online) Phase diagrams for the two sets of models generated from IUFSU (left 
panel) and TM1 (right panel). Different colors represent droplet, rod, slab, tube, bubble, and 
homogeneous phases as indicated in the legend. 
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FIG. 3: (Color online) Energy per nucleon of the pasta phase relative to that of homogeneous 
matter, A E, as a function of rib obtained from the TF (solid lines) and CP (dashed lines) calcula¬ 
tions using the original IUFSU model. Different colors represent different pasta phases as shown 
in Fig. [2j 
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FIG. 4: (Color online) Size of the Wigner-Seitz cell, r ws (thick lines), and that of the inner part, 
rd (thin lines), as a function of rib obtained in the TF approximation. The results calculated with 
the smallest L in TM1 (L = 40 MeV) and IUFSU (L = 47.2 MeV) are shown by solid lines, where 
the pasta phase structure changes from droplet (green) to rod (red), slab (blue), tube (violet), and 
bubble (cyan) as the density increases. For comparison, the results with the largest L in TM1 
(L = 110.8 MeV) and IUFSU (L = 110 MeV) are shown by green-dashed lines, where only the 
droplet configuration appears before the crust-core transition. 
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FIG. 5: (Color online) Comparison of equilibrium properties of the inner crust between the TF 
(solid lines) and CP (dashed lines) methods using the original IUFSU model. The size of the 
Wigner-Seitz cell, r ws (thick lines), and that of the inner part, ra (thin lines) (a), the electron 
chemical potential fj, e (b), and the average proton fraction Y p (c) are plotted as a function of n^. 
The different colors represent different pasta phases. 




FIG. 6: (Color online) Equilibrium properties of the Wigner-Seitz cell obtained in the TF approx¬ 
imation for modified versions of TM1 (upper panel) and IUFSU (lower panel) with several values 
of L. The average proton fraction Y p (a), the neutron densities of the liquid phase and the gas 
phase, n n ^L and n n> G, at the center or boundary of the cell (b), and the proton density of the liquid 
phase, n Pt L, at the center or boundary of the cell (c) are plotted as a function of rib- Small jumps 
are observed at the transition between different pasta shapes. 


22 






























0.00 0.02 0.04 0.06 0.08 0.10 

n b (fnV 3 ) 


0.00 0.02 0.04 0.06 0.08 0.10 


n b (fm- 3 ) 


0.00 0.02 0.04 0.06 0.08 0.10 

n b (fm- 3 ) 


FIG. 7: (Color online) Same as Fig. 0 but for chemical potentials of electrons fi e (a), neutrons fi n 
(b), and protons n p (c). 
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FIG. 8: (Color online) Density distributions of neutrons (solid lines) and protons (dashed lines) 
in the Winger-Seitz cell at rib = 0.04, 0.06, 0.07, 0.08, and 0.09 fm -3 (top to bottom) obtained in 
the TF approximation. The results with the original IUFSU model (L = 47.2 MeV) are shown by 
thick black lines. For comparison, the results with the largest L in the set of IUFSU (L = 110 
MeV) are shown by thin red lines in the top panel. The cell boundary is indicated by the hatching. 
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L (MeV) L (MeV) L (MeV) 

FIG. 9: (Color online) Crust-core transition density n^t (a), proton fraction at the transition 
point, Ypj (b), and pressure at the transition point, Pt (c), as a function of L obtained in the TF 
approximation using the two sets of models generated from TM1 (red-solid line with squares) and 
IUFSU (blue-dashed line with circles). The results obtained with the original TM1 and IUFSU 
models are indicated by the filled square and circle, respectively. 
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